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METHOD^ COMPUTER PROGRAM PRODUCT AND APPAEIATUS FOR 
ENHANCING A COMPUTERIZED TOMOGRAPHY IMAGE 

Technical Field 

The present invention relates to a methods an 
apparatus and a computer program product for improving a 
digital image, and more specifically a computerized 
5 tomography image consisting of reconstructed data, 

according to the preamble of the respective independent 
claim. 

Technical Background 

10 It is known that in computerized tomography (CT) 

images, the intensity level may be used fox 
distinguishing different types of tissue, A CT system may 
be calibrated such that the intensity values correspond 
Lo Hounsfield units (HU) as given in the table below: 





Ixitonaity 
value [HU] 


Rone ' 


1000 


Liver 


40-60 


White matter in brain 


46 


Gray matter in brain 


43 


Blood 


40 


MuscIq 


10-40 


Kidney 


30 


Cerebrospinal fluid 


15 


Water 


0 


Fat 


-50 - 100 


Air 


-1000 



The Hounsfield units normally range from -1000 to 
3000, but in some applications the HU scale is shifted 
into values between 0 and 4000. 
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As is apparent from the table above, examination of 
different types of tissue requires the inspection of 
different intensity value ranges, also referred to as 
intensity windows. 
5 Due to the limitations of display screens and also 

to limitations of the human eye, it is common to djisplay 
a maximum of 256 (2®) different gray values or shades on 
the display screen- Hence, in order to obtain a useful 
view of a CT image showing tissues having intensity 

10 values close to each other, it is necessary to display 
only a limited portion of the total HU scale. 

In the following, the displayed or processed 
intensity value range will be referred to as an 
'intensity window'. If e.g. bones are examined one 

13 typically chooses a very wide window of 1000-2000 HU, 

whereas for the examination of soft tissue, more narrow 
windows of 4 00-600 HU are used and for brain examinations 
very narrow windows of 50-100 HU are common. 

Due to the different anatomies viewed and the very 

20 different intensity ranges, the different cases require 
different noise reduction and contrast enhancement 
treatment. However it is still desirable to store the 
result in one image and be able to switch between 
intensity windows in case different tissues are to be 

25 examined. 

In US 5,594,767, an original CT image is combined 
with a smoothed version of the original CT image, the 
combination being based on a classification map- The 
smoothed image is non-selectively smoothed, i.e. it 

30 smoothes the entire image, regardless of the values of 
the respective image elements . 

Summary of the Invention 

It is an object of the present invention to provide 
35 an improved method for improving a CT image. 

The objective is wholly or partially achieved by a 
method and a computer program product according to the 



004 03/22 MON 15:46 FAX 



(21005/026 



3 

appended independent claims. Embodiments are set forth in 
the appended dependent claims and in the following 
description. 

Hence ^ there is provided a method for enhancing a CT 
5 image composed of a plurality of elements, each having an 
intensity value indicative of a tissue type. The method 
comprises receiving a first CT image, providing a second 
CT image, which substantially constitutes a copy of said 
first CT image, subjecting said second CT image to an 

10 enhancement processing, and combining said processed 
second CT image with said first CT image, whereby an 
enhanced CT image is provided, said combining being based 
on a classification with respect to intensity values of 
regions within said first and second CT images. In the 

15 method, said enhancement processing enhances said second 
CT image within a predetermined intensity value range. 

Pi CT image may be any representation of CT scan 
data, e.g. a two-, three- or four dimensional data set. 
The second CT image constitutes an enhancement 
20 processed full or partial copy of the first CT image, 
e.g- it may comprise a subset of the elements of the 
first CT image. 

The enhancement processing may be any processing 
with a view to improving the quality of the second CT 
25 image. 

According to yet another aspect, there is provided a 
method for enhancing a digital image composed of a 
plurality of elements-, each having an intensity value 
associated therewith. The method comprises receiving a 
30 first digital image, providing, by enhancement processing 
based on said first digital image, a second digital 
image, and combining said second digital image x^xth sand 
first digital image, whereby an enhanced digital image is 
provided. The combining is based on a classification with 
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respect to intensity values of regions within said first 
and second digital images. The enhancement processing is 
performed with respect to a predetermined intensity value 
range. The digital image may be any type of medical 
3 image, such as a CT image, an x-ray image etc. The 

digital image may be two-, three-, or four dimensional. 

This approach allows for intrinsically different 
processing of different intensity windows. Also, the 
above described method enables the use of different 

10 processing settings for different intensity windows and 
then combines the results into one image . 

The presented approach also allows the use of more 
advanced processing methods that for instance adapt 
locally to the image structure in an anisotropic way. 

15 According to further aspects, there are provided a 

computer program product and an apparatus for 
implementing the above described method. 

Brief des cription of the drawings 
20 Embodiments will now be described in more detail 

with reference to the appended drawings - 

Fig- 1 is a schematic flow chart illustrating a 
method according to the present disclosure- 

Fig. 2 is a schematic flow chart illustrating a 
25 detail of the flow chart of Fig. 1. 

Fig. 3 is a schematic diagram of a device 1 in which 
the method of Figs 1 and 2 may be implemented. 

Fig. 4 is a section of a simplified CT image - 
Fig- 5 is a first processed version of the CT image 
30 in Fig. 4. 

Fig. 6 is a second processed version of the CT image 
in Fig. 4 . 

Fig- 7 is a combined CT image based on the first and 
second CT images. 
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Figs 8-10 illustrates masks used in combining the CT 
images . 



Description of Embodiments 

5 For a description of the general function of a CT 

system, reference is made to US 5,594,767. Hence, the 
description of how a CT image is generated is left out. 

Fig. 1 flow chart illustrating a method for 
improving CT scan data- In step 101 , a CT image is 

10 received. It is understood that a CT image, for the 

purpose of this description, may ba received in different 
manners and formats, e.g. from a storage device, via a 
computer network, or directly from a computer that 
performed an image reconstruction based on CT detector 

15 data. 

Furthermore, in the following description, the CT 
images are assumed to be two-dimensionally indexed data 
sets, composed of a plurality of elements, each of which 
being associated with an intensity value, which may be 
20 presented in Hounsfiold units. Hence, the original CT 
image received in step 101 has been subjected to some 
known type of preprocessing and/or image reconstruction 
such as to provide a data set that is processable as an 
image . 

25 In step 102, processing parameters are received. 

Such processing parameters may include an indication of 
what intensity value range or ranges to study, an 
indication of what type or types of processing to 
perform, an indication of parameters subject to which the 

30 respective processing is to be performed. The input or 

inputs in step 102 may be received in different manners, 
e.g. in the form of a selection made by a us.er, in the 
form of an automatic determination made based on what is 
to be studied (e.g. what type of tissue), or it may be 

35 preprogrammed. 

It is recognized that steps 101 and 102 may be 
performed in the reverse order. 
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In a first embodiment, the original CT image may be 
used as input to the processing step 103, whereby an 
enhanced image is provided as an output, and whereby the 
original CT image, or a copy thereof, is maintained. 
5 Optionally, one or more CT image copies may be 

provided based on the original CT image. The CT image 
copy may be used as input to the processing step 103, so 
as to maintain the CT image intact- The CT image copies 
may be complete or partial copies of the original CT 

10 image. For example, a CT image copy may be a copy of only 
a region of the original CT image, such that the CT image 
copy becomes smaller than the original CT image. As 
another example, a CT image copy may be a copy of only 
those elements of the original CT image, x>rhich have 

15 intensity values within a certain, predetermined 

intensity value range. The CT image copies may be stored 
in a working memory or in a non- volatile memory. 

Step 102 and the step of providing the CT image copy 
may be performed in reverse order, 

20 In step 103, the original CT image and/or the CT 

image copies are subjected to enhancement processing^ 
whereby an enhancement processed CT image is obtained. 

Each one of the images may be processed in a 
predetermined manner. In one embodiment, the different 

25 images are subjected to processing within certain 

intensity windows. In another embodiment, the different 
images are subjected to different types of processing, 
i.e. to diff erent . algorithms . In yet another embodiment, 
the different images are subjected to the same type of 

30 processing but with different processing parameters. In 
yet another embodiment, one of the images is not 
subjected to processing at all, whereas the other images 
are. 

In one embodiment, the CT image is maintained and 
35 then combined with the enhancement processed CT image. In 
yet another embodiment, also the original CT image is 
subjected to enhancement processing and then combined 
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with another enhancement processed CT image. Hence/ two^ 
three or more enhancement processed CT images may be 
provided and combined - 

There are many different types of processing of the 
5 original CT image and/or the CT image copies : low-pass 
filters^ contrast enhancement using unsharp masking as 
described in US 5,038^387, rank filters such as the 
median and its variants and more advanced processing 
methods such as adaptive filtering^ such as described in 

10 Westin,C.-Fw Richolt^ J., Moharir, V. & Kikinis, R. : 
Affine Adaptive Filtering of CT Data / Medical Image 
Analysis 4, 161-177. 2000 and in US 6,556,720, Such 
adaptive filtering may for example be based on local 
structure analysis^ as described in US 4,747,150, US 

15 4,747,151, US 4,747,152 and Carmona, R.A. & Zhong, S.i 
Adaptive smoothing respecting feature directions , IEEE 
Transactions on Image Processing 7(3), 353-358- 1998. 
Other types of processing includes mean-shift filtering 
as described in Dominguez, G.F., Bischof, H. & Beichel/ 

20 R.: Fast 3P Mean Shift Filter for CT Images , SCIA, 438- 
445, 2003r G5teborg, Sweden, Springer, LNCS 2749; 
variational methods as described in Niessen, W,, Ter Haar 
Romeny, B. & Viergever, M.A. : Geodesic Deformable Models 
for Medical Image Analysis , IEEE Transaction on Medical 

25 Imaging 17(4), 634-641, 1998 and Gilboa, G., Sochen, N- fi 
Zeevi, Y,Y.: Forward and Backward Diffusion Processes for 
Adaptive Image Enhancement and Denoising , IEEE 
Transaction on Image Processing 11(7), 699-703, 2002; 
multiband techniques such as those described in Dippel, 

30 S,, Stahl, M., Wiemker, R, & Blaffert, T.: Multiscale 
Contrast Enhancement Radiographies ; Laplacian Pyramid 
Versus Fast Wavelet Transform , IEEE Transactions on 
Medical Imaging 21(4), 343-353, 2002, Vuylsteke/ P.P. & 
Schoeters, E,P. : Multiscale Image Contrast Amplification 

35 (MUSICA) , SPIE Vol, 2167 Medical Imaging: Image 

Processing, 551-560. 1994 and in US 5,644,662; or 
wavelet techniques as described in Yang, G.Z-, Rubens, 
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M.B- & Hansen, D.M. : The use of hybrid structural 
filtering for the detection of obliberative bronchiolil:is 
from computed tomographic images ^ IPA, lEE Conference 
Publication No. 443r 631-635, 1997. 
5 In one embodiment, it is determined beforehand which 

data subsets of the original CT image or CT image copy 
need to be processed, whereby all other data elements are 
excluded from the processing, in order to increase the 
overall execution speed. However, in a simple 

10 implementation the entire data set can be processed in 
parallel, or sequentially- 

There is no explicit need to segment the various 
intensity windows beforehand and it suffices to apply the 
various enhancement schemes adapted to each window in 

15 parallel. However, in order to speed up the overall 

calculation, it is possible? to compute masks (weights) 
for each window as described below and only apply the 
computation where the associated weight differs from 
zero . 

20 In step 104, the enhancement processed images and 

possibly also the original CT image (or a copy thereof) 
are combined or merged such that a combined CT image is 
obtained. Merging of the results obtained by the 
different processing methods or parameter settings may be 

25 based on the associated intensity window or intensity 
value range- As changes are more visible for smaller 
window widths/ the results may be ranked according to the 
window width, such that smaller windows are given 
priority over larger windows . 

30 In the following^ a description of an algorithm for 

merging the images will be given with reference to Fig- 
2, wherein the windows are assumed to be specified by 
their respective center c and width w. 

The algorithm may consist of the following steps, 

35 which are performed for all data elements: 

In step 1041, the CT images may be prioritised according 
to window width (intensity value range) , such that the 



r^^w^xt MfMfirlArl h\r i IQDTn frnm tho IPW imartc^ nafsihai^e on Od/05/2005 



004 03/22 MON 15:46 FAX 



@011/02S 



9 

window having the smallest width is placed first: wl < w2 
< < wN. By prioritizing the intensity windows, it is 
possible to control which window is to be used in case of 
overlapping areas. As an alternative, the prioritizing 
5 may be based on a user input, whereby the user determines 
which intensity window is to be used. 

In step 1042, regjon masks are determined by 
thresholding the original image, according to the 
following exemplifying pseudo code, wherein maskS is the 
10 thus far segmented data, orglm is the original CT image, 
maskx is a mask for the respective CT image copy and 
maskB is a mask for the non-segmented parts of the 
original image - 

masks = 0 
15 for X = 1 ... N 

{ 

maskX = 1 if (cX - wX/2 ^ orglm cX + wX/2) and (masfcS - 
0) 

masks — masks 4- maskX 
20 ) 

masicD = non-segmented parts 

Step 1042 may also be expressed according to the 

f ollovjing : 
25 1. initialize: M^(m^yO-0 \fn,m 

x-1 



[0 else 



30 4. if X < N then x = x-i-1 goto step 2 

where Mi indicates mask number i and i=X/ s, b, I 
represents the original image data and m,n indexes the 
35 data elements. 
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It is recognized that steps 1041 and 1042 may be 
performed in the reverse order . 

In step 10 44^ which is optional/ the masks are smoothened 
by convolving with a low-pass kernel designated G: 
5 masJcX » G maskX for X « 1 ... N/ B 



Step 1044 may also be expressed according to the 
following: 

10 which Is to be \inderstood as convolution of the entire 
mask array- 
In step 1045/ the region masks are normalized according 
to the following exemplifying pseudo code: 

maskX = maskX / (maskl + mask2 + ... + xnaskN + maskB) for X « 1 



Step 1045 may also be expressed as; 
Mj.(^m,n) =■ — — 



20 

In step 1046, the results are combined accorciing to the 
following exemplifying pseudo code, wherein enhim designates 
the enhanced combined CT image and wherein resultX etc 
indicates the output of the processing step 104a-c for CT 
25 image copy nxamber X: 

enhlm = maskB • orglm + maskl - result 1 + ... + maskN • resultN 

Step 1046 may also be expressed as: 
3 0 EOfh w) = -A/tf (w, w) • I{m, ri) + M, (iw, w) • Jf, (w, «) + ... + Af (/w, w) ■ jR^ ( w, «) 



where E denotes the final enhanced image and Rir i=^l-.Nr 
the previously obtained results from the processing 
according to the individual intensity windows- 
35 It should be noted that in order to avoid artifacts 

near the boundaries of the respective regions, it is 
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possible to apply a low-pass filter to the segmented 
masks in step 1044 above. In one embodiment, a small 
Binomial or Gaussian filter may be chosen/ however any 
other low-pass kernel will do, 
5 The method segments the original data and then 

applies a spatial transition (smoothing) between areas to 
smooth the transition. 

The combining as described above ^ and in particular 
the normalizing^ may be useful as the proposed method 

10 recombines intensity values and it is essential that the 
local mean of these Hounsfield units is not altered by 
the processing. 

In step 105, the combined CT image ts output/ e.g. 
by being displayed on a screen, stored in a memory or 

15 sent via a network- 
Optionally, a step 1047 comprising a morphological 
closing and/or opening may be included with respect to 
the masks r thereby removing or filling small gaps in the 
respective masks. The morphological closing and/or 

20 opening step 1047 may be arranged to take place at 

anytime between the determination of the region masks 
1042 and the combining 104 6 of the CT images - 

Fig. 3 is a schematic diagram of a device 1 in which 
the above described method may be implemented. The device 

25 1 comprises receiving means 2 for receiving a CT image 
from e.g. a network, a storage medium or from a 
preprocessing device. The device 1 further comprises a 
processing unit 3, which is adapted for performing the 
above described method or methods. The device 1 also 

30 comprises output means 4 for outputting an enhanced CT 

image. The processing unit may consist of a programmable 
processor/ which by means of software is programmed to 
perform the described method. Tn another embodiment, the 
preprocessing device may take the shape of a digital 

35 signal processor that xs adapted to perform the described 
method. Such a preprocessing device may be provided in 
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the form of e.g. an ASIC or any structure having similar 
functionality - 

It is also possible to implement the method in the 
form of a propagated signal comprising components for 
5 performing the steps of the above described method when 
said components are executed in a computer - 

Below, the processing of a section of a simplified 
2D CT image will be described with reference to Figs 4- 
10. In Figs A-1 r the horizontal axis illustrates image 

10 elements in e-g. a vertical line of the respective CT 

image. The vertical axis illustrates the HU value of the 
respective image element in the section. 

Fig. 4 is a section of an original CT image^ e-g- as 
received in step 101 {Fig. 1) . In the described example, 

15 two intensity windows will be considered: a wide window 
ranging from about 1250 HU to about 3650 HU, which is 
illustrated by the dotted lines in Fig. 4; and a narrow 
window ranging from about 2150 to about 2600 HU, which is 
illustrated by the dashed lines in Fig- 4. The respective 

20 window may be determined as described above. In one 

particular embodiment, the limits of the intensity window 
may be determined by a graphical user interface, wherein 
the lines illustrating the limits of the intensity window 
are presented in relation to a section of the image that 

25 is to be processed, whereby the user, by means of a 

pointing device, drags the lines to desired positions. 

Fig- 5 is a first processed version of the original 

CT image in Fig. 4, as provided by step 103 (Fig. 1) , In 
the illustrated embodiment, the original CT image of Fig. 
30 4 has been subjected to a low-pass filter, with respect 
to the wide window. Such low-pass filtering suppresses 
noise - 

Fig. 6 is a second processed version of the original 
CT image in Fig. 4. In the illustrated embodiment, the 
35 original CT image ties been subjected to a structure- 
adaptive processing with respect to the narrow window. 
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Such structure-adaptive processing preserves small 
details in the CT image, but does not suppress noise to 
the same extent. 

Fig. 7 is a combined CT image based on the first and 
5 second processed CT images, as provided by step 104 (Fig. 
1) . The combination has been made with respect to the 
intensity windows associated with the first and second 
processed CT images. 

Figs 8-10 illustrates masks used in combining the CT 

10 images* Fig. 8 illustrates the mask associated with the 

narrow window, i.e. with all image elements having values 
within that window. Fig- 9 illustrates the mask 
associated with the wide window, i-e, with all image 
elements falling within that window, but outside the 

15 narrow window. Fig- 10 illustrates the background mask, 
i.e. with all image elements falling outside both the 
narrow and the wide window. It should be noted that the 
masks have a slightly smooth transition, due to the low- 
pass filtering of the original binary masks. 

20 It should also be noted that in case a morphological 

closing and/or opening step 1047 (Fig. 2) would be 
performed with respect to the masks illustrated in Figs 8 
and 9, then the unevenness shown in the regions between 
50-90 and 170-210 on the horizontal scale would be 

25 removed. 

Referring back to Fig- 7, it is noted that for those 
elements whose HU values fall outside both windows, the 
data from the original CT image (Fig. 4) is used. For 
those elements whose HU values fall v/ithin the narrow 
30 window, the data from the second CT image (Fig. 6) is 

used. For those elements whose HU values fall within the 
wide window, but outside the narrow window^ the data from 
the first CT image (Fig- 5} is used- 

•'M *• • II ,••!,;!•:/:/:• *cr 

. I."- 
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Even though the above description focuses on 2D 
image slices^ the method is analogously applicable to the 
processing of higher dimensional data such as direct 3D 
or 4D enhancement of CT data. 
5 It should also be noted that although the 

description focuses on CT images^ the described method is 
analogously applicable to any type of image, and in 
particular to medical images ^ such as x-ray^ MRI, etc. 

Additional areas of application for the method 

10 described above include, but is not limited to, non- 
destructive testing, e.g. for use in quality assurance in 
the manufacturing industry; analysis of logs in saw mills 
and in connection with airport luggage security checks. 



15 
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CLAIMS 

1. A method for enhancing a CT image composed of a 
plurality of elements^ each having an intensity value 
indicative of a tissue type^ the method comprising: 

receiving (101) a first CT image, 

providing, by enhancement processing (103) based on 
said first CT image, a second CT image, and 

combining (104) said second CT image with said first 
CT image/ whereby an enhanced CT image is provided, 

said combining being based on a classification with 
respect to intensity values of regions within said first 
and second CT images, 

characterized in that 

said enhancement processing (103) is performed with 
respect to a predetermined intensity value range. 

2. The method of claim 1, further comprising 
receiving an indication of said predetermined value range 

20 and associating said predetermined intensity value range 
with said second CT image - 

3. The method of claim 1, wherein said enhancement 
processing (103) is adaptive to a local structure defined 

25 by at least some of said plurality of elements. 

4. The method of claim 3, wherein said local 
structure is defined by a group of elements whose 
intensity values are within said predetermined intensity 

30 value range. 
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5. The method of claim Ir wherein said enhancement 
processing (103) comprises applying a non-linear filter 
to said second CT image . 

5 6. The method of claim 1, wherein said enhancement 

processing (103) is selected from a group consisting of a 
noise reduction using a low pass filter^ a contrast 
enhancement using unsharp masking, a rank filtering^ an 
adaptive filtering, a mean-shift filtering, a variational 
10 method, a multiband technique and a wavelet technique - 

1 . The method as claimed in any one of the preceding 
claims, wherein combining (104) said first and second 
digital CT images comprises: 

15 determining (1042) a first region mask for said 

first CT image^ said first region mask defining an area 
within the first CT image, whose elements have intensity 
values within a first intensity value range, 

determining (1042) a second region mask for said 

20 second CT image, said second region mask defining an area 
within the second CT image, whose elements have intensity 
values within said predetermined intensity value range, 
and 

combining (1046) said first and second CT images, 
25 weighted by their respective region masks, whereby said 
enhanced CT image is provided, 

8. The method of claim 7, further comprising 
prioritizing (1041) said first and second CT images, 
30 whereby an element of a CT image having a higher priority 
is included in the enhanced CT image and a 
correspondingly located element of a CT image having a 
lower priority is excluded from the enhanced CT image . 

. V:* ."'M.:". • -A* 
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9. The method of any one of claims 7 or 8r further 
comprising smoothing (1044) said first and second region 
masks . 

5 

10. The method of any one of claims 7-9, further 
comprising normalising (1045) said region masks, 

11. The method of any one of claims 7-10, further 
10 comprising subjecting at least one of said first and 

second region masks to a morphological closing and/or 
opening algorithm. 

12. The method of any one of the preceding claims, 
15 wherein said first CT image is selected from a group 

consisting of a two-dimensional array, a three- 
dimensional array and a four-dimensional array. 

13- The method as claimed in any one of the 
20 preceding claims, wherein said first CT image is 

subjected to a second enhancement processing prior to 
said combining (104) . 

14. The method as claimed in claim 13, wherein said 
25 second enhancement processing is performed with respect 
to a second predetermined intensity value range. 

35. A computer program product comprising software 
code portions for performing the steps of any one of 
30 claims 1-14, when said product is run on a computer - 

16. A storage medium having stored thereon a 
computer program product according to claim 15. 
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17. A propagated signal comprising components for 
performing the steps of any one of claims 1-14. 

18. A device for enhancing a CT image composed of a 
5 plurality of elements , each having an intensity value 

indicative of a tissue type, the device comprising 

receiving means (2) for receiving a first CT image ^ 
processing means (3) arranged for providing, by 

enhancement processing (103) based on said first CT 
10 image, a second CT image, and 

for combining (104) said second CT image with said 

first CT image, whereby an enhanced CT image is provided, 
said combining being based on a classification with 

respect to intensity values of regions within said first 
15 and second CT images, 

characterized in that 

said processing means (3) is adapted for enhancement 
processing with respect to a predetermined intensity 
value range. 

20 

19. A method for enhancing a digital image composed 
of a plurality of elements, each having an intensity 
value, the method comprising: 

receiving (101) a first digital image, 
25 providing, by enhancement processing (103) based on 

said first digital image, a second digital image, and 

combining (104) said second digital image with said 
first digital image, whereby an enhanced digital image is 
provided, 

30 said combining being based on a classification with 

respect to intensity values of regions within said first 
and second digital images, 

characterized in that 
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said enhancement processing (103) is performed with 
respect to a predetermined intensity value range. 
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ABSTRACT OF THE DISCLOSURE 

A method is disclosed for enhancing a CT image 
composed of a plurality of elements ^ each having an 
5 intensity value indicative of a tissue type. The method 
comprises receiving (101) a first CT image, providing, by 
enhancement processing (103) based on said first CT 
image, a second CT image, and combining (104) said second 
CT image with said first CT image, whereby an enhanced CT 

10 image is provided. The combining is based on a 

classification with respect to intensity values of 
regions within said first and second CT images. The 
enhancement processing (103) is performed with respect to 
a predetermined intensity value range. 

15 The method may be implemented as a computer program 

product or as a device, 
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